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Compared with classical search algorithms, Grover quantum algorithm [Phys. Rev. 
Lett., 79, 325(1997)] achieves quadratic speedup and Brüschweiler hybrid quantum 
algorithm [Phys. Rev. Lett., 85, 4815(2000)] achieves an exponential speedup. In this 
paper, we report the experimental realization of the Brüschweiler algorithm in a 3- 
qubit NMR ensemble system. The pulse sequences are used for the algorithms and the 
measurement method used here is improved on that used by Brüschweiler, namely, instead 
of quantitatively measuring the spin projection of the ancilla bit, we utilize the shape of 
the ancilla bit spectrum. By simply judging the downwardness or upwardness of the 
corresponding peaks in an ancilla bit spectrum, the bit value of the marked state can 
' be read out, especially, the geomètric nature of this read-out can make the results more 

I ' robust against errors. 

PACS numbers: 03.67.Lx, 82.56.Jn, 76.60.-k. 
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I. INTRODUCTION 



' Quantum algorithms are very important in quantum computing. One can find this point in Deutsch 
y—\ , and Josza's quantum algorithm which demonstrates the incomparable advantage of quantum computing 
|EJ. Two more famous quantum algorithms which are closely related to practical applications of quan- 
tum computation are: Shor's factoring algorithm j| and Grover's quantum search algorithm j^]. The 
■ factorization of a large number into prime factors is a difficult mathematical problem because existing 
f || classical algorithms require exponential times to complete the factorization in terms of the input. How- 
! 1 > . ever, Shor's quantum algorithm drastically decreases this to polynomial times. Another similar example 
is searching marked items from an unsorted database. Actually, many scientific and practical problems 
can be abstracted to such search problem. Hence it is a very important subject. Classically, it can only 
be done by exhaustive searching. Unlike Shor's algorithm, Grover's quantum algorithm achieves only 
quadratic speedup over classical algorithms, namely the number of searching is reduced from O(N) to 
0(\f~N). However, it has been proven that Grover's algorithm is optimal for quantum computing B. 
The strong restriction of the optimality theorem can be broken off if we go out of quantum computation 
and then exponential speedup may be achieved. Using nonlinear quantum mechanics, Abram and Lloyd 
H have constructed a quantum algorithm that achieves exponential speedup. However the applicability 
of nonlinear quantum mechanics is still under investigation, let alone the realization of their algorithm 
at present. 

Recently, by using multiple-quantum operator àlgebra, Brüschweiler put forward a hybrid quantum 
search algorithm that combines DNA computing idea with the quantum computing idea |6|] . The new al- 
gorithm achieves an exponential speedup in searching an item from an unsorted database. It requires the 
same amount of resources as effective pure state quantum computing. There are several known schemes 
for quantum computers, such as cooled ions 0, cavity QED ||, nuclear magnètic resonance |9j and so on. 
NMR tcchniquc is sophisticated and many quantum algorithms have been realized by using NMR system 
fï(j| |Ï6| . Many studies show that the NMR system is particularly suitable for the realization of such algo- 
rithms, in which ensembles of quantum nuclear spin system are involved. Strictly speaking, Brüschweiler 
algorithm is not a pure quantum algorithm, and thus the realization of Brüschweiler algorithm in NMR 
enjoys the freedom from the debaté P^íq] about the quantum nature of the NMR computation using 
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efïective pure state. Because this algorithm is exponentially fast, it tales much shorter time to finish a 
search problem and this also makes the algorithm more robust again decoherence. 

In this paper, we report the experimental realization of Brüschweiler algorithm in a 3-qubit homo- 
nuclear systcm. In the procedure, we have improved the measurement method used by Brüschweiler in 
his paper ||. Instead of measuring the ancilla bit's spin polarization, we utilize the shapes of ancilla bit's 
spectra, i.e. by judging the downwardness or upwardness of the corresponding peaks in the spectrum, 
the bit value of the marked state can be read out. Since the geomètric property of the spectrum is easy 
to be recognized, this makes the algorithm more tolerant to errors. Our paper is organized as follows. 
After this introduction, we briefly describe Brüschweiler's original algorith m i n section iïïïand then we 
introduce our modification part based on Brüschweiler algorithm in section III. In sectionjiV|, we present 
the details of the pulse sequences of the algorithm and the results of our experiment. Finally, a summary 
is given. 



II. BRÜSCHWEILER'S ALGORITHM 



NMR techniques lie far ahead of other suggested quantum computing technologies. However, during 
recent years, the ràpid developing tendency becomes slow and slow. People have taken more effort in 
preparing an effective pure state, but compared with a pure state quantum computer, there is no essential 
speedup. Brüschweiler's wonderful idea may shed a light on this area, he takes advantage of the mixed 
state nature in the NMR system and achieves an exponential speedup in searching an unsorted database. 
For convenience in following discussion, here we repeat the main idea of the Brüschweiler's algorithm in 
brief (in detail, see Refs. §J2ÏJ). 

As is well-known, the preparation of the effective pure state is one of the most troublesome part in a 
NMR quantum computing experiment. On the other hand, the effective pure state also sets a restriction 
on the number of qubits [|19|j20| . The effective pure state is represented by the density operator 



p = (1 - e)2~ n i + e |00...0) (00...0I . 
At room temperature, under the high temperature approximation we have 

nhv 
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In Eq. ([!]), the second term's contribution to the outcome is scaled by the factor e, which decreases 
exponentially with n, namely, the number of qubit, but the first term has no contribution at all |22[ . 

In NMR ensemble system, the state can be represented by density operators which are linear combi- 
nations of direct products of spin polarization operators P]2l|1. In a strong external magnètic field, the 
eigenstates of the Zeeman Hamiltonian 



\<f> in ) = 1001—01) = \aa/3...ap) , 
are mapped on states in the spin Liouville space 
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represent respectively spin up and spin down state of the spin. Usually, the oracle or query is a computable 
function /: f(x) = for all x except for x = z which is the item that we want to find out for which 
f(z) — 1. Usually, the oracle can be expressed as a permutation operation which is a unitary operation 
Uf, implemented using lògic gates ||. In Brüschweiler algorithm, an extra bit(also called the ancilla bit) 
is used and its state is represented by Iq. The output of the oracle is stored on the ancilla bit I whose 
state is prepared in the a state at the beginning. The output of / can be represented by an expectation 
value of Iqz for a pure state 
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/ = F(I^a in ) = i - Tr(U f I^a m U+I 0z ). (7) 

If <7j„ happens to satisfy the oracle, then Iq is changed to Iq. This gives the value of the trace equal 
to —1/2, and hence / equals to 1. The input of / can be an mixed state of the form p = XX'=i a ó 
where <jj is one of the form in Eq. (^): 

- N N-ï 

f = Ew ^) = nE J o ^) + (8) 

The oracle is applied simultaneously to all the components in the NMR ensemblc. The oracle operation 
is quantum mechanical. Brüschweiler put forward two versions of search algorithm. Wc adopt his second 
version. The essential of the Brüschweiler algorithm is as follow: suppose that the unsorted database has 
TV = 2™ number of items. We need n qubit system to represents these 2 n items. The algorithm contains 
n oracle queries each followed by an measurement: 

(1) Each time, IqI% (k = 1,2, n) is prepared. In fact, the input state Iq...1...I%...1... is a highly mixed 
state 21 1. In the following text, the identity operator will be omitted. This Liouville operator actually 



represents the 2™ 1 number of items encoded in mixed state: 



E 
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= E 1*1*2 — Ík-lOÍk+1— i n ) {hh— ik-iQik+i— i n \ ■ (9) 

Í1,Í3,— ,ín=0,l 

This mixed state contains half of the whole items in the database. The fc-th bit is set to a. The othcr 
half of the database with fc-th bit equals to (3 (or 1) is not included. (2) Applying the oracle function to 
the system. As seen in eq. (Q), the operation is done simultaneously to all the basis states. If fc-th bit 
of the marked state is 0, then the marked state is contained in eq. (g). One of the 2 n terms in equation 
(^) satisfies the oracle and the oracle changes the sign of the ancilla bit from a to 0. If one measures the 
spin of ancilla spin after the function /, the value will be / = (2™ — 1) x 1/2 + 1/2 — (2" — 2) x 1/2 = 1. 
If the fc-th bit of the marked state is 1, then the state @ will not contain the marked item. Upon the 
operation of the function /, there is no fiip in the ancilla bit. A measurement on the ancilla bit's spin Iq z 
will yield / = 1/2 x (2" — 1) + 1/2 — (2") x 1/2 = 0. However, without obtaining the value of /, we can 
know the marked state by measuring the ancilla bit's spin. If one measures the spin of ancilla spin after 
the oracle, the value will be (2 n_1 - 1) x 1/2-1/2 = A/4 — 1 for the fc-th bit of the marked state being 
0. If the fc-th bit of the marked state is 1, then the state (|9|) will not contain the marked item. Upon the 
operation of the oracle, there is no fiip in the ancilla bit. A measurement on the ancilla bit's spin Iq z will 
yield 1/2 x (2™ _1 ) = A/4. Therefore by measuring the ancilla bit's spin, one actually reads out the fc-th 
bit of the marked state. (3) By repeating the above procedure for k from 1 to n, one can find out each 
bit value of the marked state. 

In the following, we give a simple example with N = 4 for illustrating the algorithm, and the cxample is 
realized in an experiment. The example is used for demonstration. The advantage of the algorithm will be 
seen if the number of qubit becomes large. Suppose the unsorted database with four items {00, 01, 10, 11} 
is represented by Zeeman eigenstates of the two spins I±, 12- The item z = 10 is the one which we want. 
That is to say, / = 1 for z = 10, which is expressed as 1x1% ■ For the other three items, {00(1" 1%), 
01(/f/|), ll(lfif)}, / = 0. Function / can be realized by a permutation illustrated in Fig.l (similar as 
Figure 2 in Ref. pj). The extra qubit Iq is included in the permutation. 

First we prepare a mixed state IqIi, which is the sum of I0I1I2 + IqIiIÍ- Then the permutation 
described in Fig.l is operated on this mixed state. Since the first bit of the marked state is 1, the 
permutation will have no effect on the ancilla bit because it is obviously that the state IqI" will not 
contain the marked state. I^Ifl^, IqI" 1% each contributes 1/2 to the spin of the ancilla bit. Upon 
measurement of the ancilla bit on its spin, the intensity will be 2x1/2 = 1 unit. That telis us that the first 
bit of the marked item is 1( in state Jf ). Secondly, we prepare another state, IqI^ — Iq!?^ + ^0-^1^2- 
We get output I^I^Ig + I^IjI^ after the action of permutation /. Measuring the spin of ancilla bit, 
we get 0, since IqI^I^ and I^I^I^ contribute to the spin measurement equally but with opposite signs. 
Then this telis us that the second bit is 0( in state Jf). After these two measurement, we have obtained 
the marked state. In the actual experiment, we have modified the measuring part of the algorithm. We 
read out the bit vàlues by looking at the shape of the ancilla bit. It is more clearly and concise. 
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III. MODIFICATION TO THE ORIGINAL ALGORITHM 



Need not measure the I§, we can distinguish the state of the ancilla bit by the shape of its spectrum. 
Because difïerent initial state InK has the same form, except for the difference in the k subscript, it 



o will have similar shapes for Iq 1% and Iq Ifc . We use the shape of the 



is natural that the spectrum /( 

spectrum of the state Iq 1% as a reference where k = 1, 2 • • •. First, the phase of Iq Jf is determined as 
making peaks of the spectrum up. In this NMR system, the lo bit has J coupling to both I\ and I2 and 
there are only two peaks in the Iq spectrum for the state IqI", the localities of two peaks are determined 
by the order of the nuclei, which is the first, second- • • The lo spectrum of IqI" before the operation of 
the permutation Uf is given in Fig. 2(a). After the permutation operation, we measure the spectrum of 
Iq in new states again. If the shape of the spectrum is the same as the one before the oracle, i.e., two 
peaks are up still, then the permutation operation has not changed the state IqI^i and this means that 
Ik is 1, that is to say, the fc-th bit value of the marked states z is 1. If the fc-th bit of the marked state 
is 0, the ancilla bit will flip after the operation of the permutation Uf. We can see from density matrices 
before and after the query operation Uf. Before the query is evaluated on the mixed state IqI", the 
density matrix(apart from a múltiple of the identity matrix and a sealing factor) is 
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When we measure the spectrum of ancilla bit I , the left peak, corresponding matrix element 51 and the 
right peak, corresponding the matrix element 62, do not change. This indicates that the shape of the 
spectrum does not change. As for the second step, before the query is evaluated on the mixed state IqI^, 
the outeome matrix is 
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and after the query, the matrix becomes 
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The left peak( (51) matrix element) does not change, but the right peak, ((72) matrix element) changes 
sign. Thus the right peak of the spectrum will be downward. 

This method of "reading out" the bit of the marked state is effective. Since it depends on the shape of 
the spectrum, a topological quantity, it is insensitive to errors as compared to the quantitative measure- 
ment of the spin of the ancilla bit. 



IV. THE REALIZATION OF THE ALGORITHM IN NMR EXPERIMENT 



We implemented Brühweiler algorithm in a 3 qubit homomiclear NMR system. The physical system 
used in the experiment is 13 C labeled alanine li C x H 3 - 13 C°H(NH+) - 13 C 2 OOH. The solvent is D 2 0. 
The experiment is performed in a Bruker Avance DRX500 spectrometer. The parameters of the sample 
were determined by experiment to be: J02 = 54.2Hz, Jqi — 35.1Hz and J12 = 1.7Hz. In the experiment, 
1 H is decoupled throughout the whole process. 13 C°, 13 C 1 and 13 C 2 are used as the 3 qubits, whose 
state are represented by Iq, li, h respectively. 13 C° is used as the ancilla bit and the result of the oracle 
is stored on it, 13 C 1 and 13 C 2 are the second and third qubit respectively. We assume the marked item 
is 10. 

Firstly, the state I^If is prepared. It is achieved by a sequence of selective and non-selective pulses, 
and J-coupling evolution. We begin our experiment from thermal equilibrium state. This thermal state 
is expressed as, 

a(0_)=/ z ° + / z 1 +4 2 . (14) 

The input state Iq!" can be written as 5(^1 + I® + l\ + 21® 1^), The identity operator does not 
contribute signals in NMR, and a scale factor is irrelevant, thus IqI" is equivalent to I® + l\ + 21® l\. 
The pulse sequence plpSpBI 



(Í) v ^ Grad ^ (4). ^ T ^(e)- y ^ Grad (15) 



applied to the thermal state produces this input state 

a(0 + ) = ^(I a z +ll+2I°X). (16) 

However in our experiment only the spectrum of Iq is needed, and only J coupling between qubit and 1 
is retained, a simplified pulse sequence is actually used in the present experiment to prepare an equivalent 
input state: 





2 ) y ' V 2 ) : V 4 



^ %m ^Grad. (17) 



Here the subscripts denote the directions of the radio frequency pulse, and the superscripts denote the 
nuclei on which the radio frequency are operated. Two numbers at the superscript mean that the pulse are 
applied simultaneously to two nuclei(In actual experiment, the pulses are applied in sequence. Because 
the duration of the pulse is very short, they can be regarded as simultaneous) . Grad refers to applying 
gradient field. r = 1/(2 Jqi) or t' = 1/(4 Jqi) is the free evolution time during which nuclear 13 C 2 is 
decoupled. The second pulse sequence is operated more easily, because only selective to Jo is considered. 
Pulse sequence ( |Ï7| ) transforms the thermal state ( |Ï4| ) into 

<Ko+) = + il + zipl) + \il + il (18) 

States (^6|) and ( |Ï8| ) are equivalent, because |/J and Z 2 does not contribute to the 1° spectrum, and the 
scaling factor does not matter. 

The oracle, represented as a permutation / is applied to this initial state: Iq I". Thcn result of the 
oracle operation is stored on the ancilla bit lo, that is, the state of the 13 C° indicates the state of the 
first bit of the marked item. Specifically the expression of the unitary operation corresponding to the 
permutation / is 
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The permutation U / can be completed using a sequence lògic gates given in Fig. 3. The left one is 
CNOT gate and the right one is Toffoli gate. The pulse sequence can be found out in Ref. ||^,|2^]. The 
pulse sequence will be very complex if we write according to the network although it is very rigorous. 
Since we assume that there is only one marked state and only the spectrum of Iq is needed, the function 
of the Uf can be realized by the pulse sequence shown below 
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where r = 1/(2 Jo i)- After the operation of the oracle, we measure the spectrum of the ancilla bit. This 
pulse sequence achieves the same result as that for the gate shown in Fig. 3: 
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Secondly, the initial state IqI^ is prepared. There are two ways to prepare this initial state. One 



method is to use a pulse sequence as in Eq. ( |ï"5| ) or 
Anothcr method is to use the swap operator in Ref. 



17|) by exchanging 1 with 2 in the superseripts. 
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onto the initial state IqI" and the state IqI^ will be obtained. In the experiment, we adopt the second 
approach. The swap operator is important in generalizing the experiment into more qubit system and 
we will discuss this later. Then we apply the permutation Uf again, and the result of the oracle is stored 
in the ancilla bit 13 C°. 

The spectra for I after the oracle query Uf operated on Iq 7f and IqI£ are given in Fig. 2 (b) and 
Fig. 2 (c) respectively. We can see clearly that the one has the same shape as the reference spectrum and 
the other one has flipped the right peak. This telis us that the first bit and the second bit of the marked 
state are 1 and respectively. Thus the marked state is 10. We also notice that there are small differences 
between the spectra before and after the permutation operations for Ifilf. These are expected due to 
imperfections caused by the inhomogeneous held, the errors in the selective pulse and in the evolution of 
chemical shift. 



V. SUMMARY 



In summary, we have successfully demonstrated the Brüschweiler algorithm in a 3-qubit homo-nuclear 
NMR system. Pulse sequences are given. A new method for "reading out" the bit value of the marked 
state is proposed and realized. The number of iteration required for this algorithm is very small. This is 
particularly propitious to resist decoherence, especially for NMR system at room temperature. Another 
advantage of this algorithm is its robustness against errors, i.e. the shape of the spectrum in reading 
the bit of the marked state has a special feature and one can easily distinguish it from others. Another 
advantage is its high probability in finding the marked state, it is 100%!. 

It should be point out that there are still several issues to be addressed in generalizing the searching 
machine to more qubit system. First, one must find a suitable molecule to act as the quantum computer. 
According to Brüschweiler's original algorithm, the ancilla qubit Iq must interact with every other qubit. 
However, in a molecule, the interaction between remote nuclear spins is very weak. This may be overcome 
by the swap operation as given in Ref. pi) . Using swap operation, we can prepare any initial state IqI% 
without the direct interaction between spin Iq and spin J^. And, the qubit also can be read out easily 
from the shape of the spectrum. All these are under consideration in our future work. 
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Q. 

FIG. 1. Representation of the oracle Uf spanning on spins li, I2 whose function corresponds to one permutation 
and results of the query stored on lo. 
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FIG. 2. Experimental realization of the Brüschweiler's algorithm using the three spins 1/2 of 13 C , 13 C-labled 
alanine dissolved in D2O. (a) The spectrum of the ancilla qubit in state before the oracle Uf, two peaks are 
upword. The shape of the spectrum is used as criterion. (b) The spectrum of the ancilla qubit after the oracle 
Uf. The first bit is 1. (c) The spectrum of the ancilla qubit in the state after the oracle Uf. The second bit is 0. 
The marked state is 10. 
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FIG. 3. Network for realizing the oracle [//. 
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